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Abstract: We have investigated the pressure dependence of the detection sensitivity of 
CO2, N 2 0 and CH 4 using wavelength modulation spectroscopy (WMS) with distributed 
feed-back diode lasers in the near infrared region. The spectral line shapes and the 
background noise of the second harmonics (2f) detection of the WMS were analyzed 
theoretically. We determined the optimum pressure conditions in the detection of CO2, 
N2O and CH4, by taking into consideration the background noise in the WMS. At the 
optimum total pressure for the detection of CO2, N2O and CH4, the limits of detection in 
the present system were determined. 

Keywords: gas sensing; wavelength modulation spectroscopy; near-infrared 
absorption spectroscopy 
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1. Introduction 

Atmospheric pollution is a problem as global warming is caused by the man-made addition of extra 
amounts of greenhouse gases such as CO2, N 2 0 and CH 4 . Current average mixing ratios of the trace 
components CO2, N 2 0 and CH 4 in the atmosphere are ~380, ~1.8 and ~0.32 parts per million by 
volume (ppmv), respectively. CO2 is released into the atmosphere when fossil fuels, coal, solid waste 
and wood products are burned, during cement production, and also when land surface cover is changed 
by humans. Increased CH 4 levels come from fossil fuels, rice cultivation, animal husbandry, biomass 
burning and landfills. The main anthropogenic sources of N2O are agriculture, and industrial sources 
including adipic and nitric acid production. Combustion of solid waste and fossil fuels also contributes 
to atmospheric N2O. The need to detect these trace gases has become increasingly important in recent 
years, both for controlling industrial processes and for monitoring air quality. Since many trace gases 
have a significant impact on the environment, the development of techniques for its fast, accurate and 
sensitive detection is required [1-5]. 

Laser based spectroscopy provides interesting advantages related to its high selectivity and 
sensitivity in the detection of trace gases. Semiconductor diode lasers emitted the near infrared (NIR) 
light have played a central role, being tunable spectroscopic light sources which exhibit a relatively low 
amplitude noise [6]. In addition, their small size and the relatively low cost have made them 
particularly suitable for the realization of transportable spectrometers employed for in-situ 
measurements [7,8]. In the NIR region, relatively weak overtone or combination vibrational transitions 
occur for trace gases relevant to the atmospheric environment. Although line strengths of these 
transitions are a few orders of magnitude lower than transitions belonging to fundamental vibrational 
bands, cavity ring down spectroscopy (CRDS) [9] and frequency modulation spectroscopy (FMS) with 
a single-mode diode laser in the near infrared region have allowed to achieve the detection of low 
concentrations (< sub ppmv) of trace gases [10]. The use of multi-pass absorption cells [11] further 
improves the sensitivity, since the effective optical path is increased up to several tens or hundreds of 
meters. Phase-sensitive techniques, which are the core of the FM techniques, significantly reduces 
the \lf electronic noise, achieving high detection sensitivity. In most cases the detection sensitivity is 
only limited by noise introduced by undesired optical fringes and drifting of laser power and detector 
sensitivity. Wavelength modulation spectroscopy (WMS) [12-16] is a kind of FMS with a modulation 
frequency lower than the spectral line width of interest. In the case of WMS, the lower modulation 
frequency allows one to use low-frequency circuits and photo detectors, reducing the complexity and 
the cost of measurement system. Recently, a high sensitive open path CH4 analyzer based on WMS 
detection and Herriott cell design with a 30 m total optical path length has been developed which is 
commercially available [17]. 

In the present study, we demonstrate the detection of greenhouse gases such as CO2, CH 4 and N 2 0 
using WMS. The aim of this work is characterized the optimum condition of the detection of 
greenhouse gases by NIR- WMS. The optimum pressure for the detection of CO2, N2O and CH4 by 
WMS is investigated both experimentally and theoretically by taking into consideration background 
noise in the WMS. 
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2. Experimental 

We used WMS in the near-infrared region for trace gas detection. Figure 1 shows a scheme of our 
experimental apparatus. All optical elements were properly designed to reduce their dimensions and 
mechanical instability as much as possible. Three distributed feed-back (DFB) diode lasers 
(NTT Electronics) with a tuning range of ±1 nm were used as light sources. The first DFB laser with a 
center wavelength of 1,572 nm was used to detect CO2 at the 3vi + V3 combination band. The second 
DFB laser with a center wavelength of 1,515 nm was used to detect the N2O line at the 3v3 overtone 
band. The last DFB laser with a center wavelength of 1,651 nm was used to detect CH4 at the 2V3 
overtone band. The laser power varied as a function of injection current with a maximum of 20 mW. 
The wavelengths of these DFB lasers were chosen for strong absorbing transitions and no interference 
of water. 



Figure 1. Schematic diagram of the experimental setup (MFC: mass flow controller. 
Func. Gen.: function generator). 
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The laser wavelengths were tuned by varying the laser temperatures, while fine tuning was 
accomplished by changing the laser diode injection currents. To minimize fringe noise from optics, 
anti-reflection (AR) coated lenses and AR coated and wedged windows were used. The output beam 
was focused by an AR coated lens (f = 50 cm) to the centre of a Herriott-type multi-pass cell [11] with 
a distance of 40.4 cm between two mirrors. In the multi-pass cell, the laser beam undergoes multiple 
reflections between two mirrors; the total number of 74 passes corresponds to the optical path length of 
29.91 m. The inner volume of the cell is 900 cm . The beam passed through the multi-pass cell was 
finally focused by a short focal AR coated lens (f = 5 cm) onto an InGaAs photodiode detector 
(Hamamatsu G5 852-1 1). 

The laser wavelength was sinusoidally modulated at 10 kHz, and scanned at 1 Hz around the 
absorption line. The laser wavelength scan and modulation were performed through a custom-made 
laser driver. The background bias signal of the laser driver used in this study was same level as that of 
a commercial driver (e.g., ILX LDC3724C). Under typical conditions, the modulated absorption signal 
was phase-sensitive detection at twice the modulation frequency (2f) using a lock-in amplifier 



Sensors 2010, 10 



4689 



(Stanford Research System SR810 DSP) with the time constant set to 3 ms. The data was acquired via 
a 16 bit AD PCMCIA card (CONTEC ADA16-8/2(CB)L) to a laptop computer and analyzed using 
software written in Lab VIEW. 

N 2 (99.99%) and premixed gases of synthesized air (79% N 2 , 21% 0 2 ), 1% C0 2 in N 2 , 1% N 2 0 in 
N 2 , and 0.1% CH 4 in air were purchased from Taiyo Nippon Sanso Co. The flow rate of the feed was 
controlled with mass flow controllers (Kofloc MODEL3660). The total pressure was monitored with a 
capacitance manometer (Setra Systems MODEL720). All measurements were performed at room 
temperature (293 ± 3 K). 

3. WMS Theory 



The general theory of the FMS and WMS has been well summarized in the works of previous 
researchers [13-15]. In the case of sinusoidal modulation with modulation index /? and frequency co m , 
the optical field E with carrier frequency a>o is given by: 



E(t) = E 0 exp[i(cy + /?sin a>j)] = E 0 exp(ico 0 t) £ /„ (j3)exp(ina> m t) 



(1) 



where the expansion in a series of nth-order Bessel functions J n characterizes the frequency 
components of the modulated light spectrum. FMS is classified by its modulation frequency co m . WMS 
uses co m that is much less than the spectral line width r of the absorption line of interest. Conversely, 
narrowly-defined FMS uses co m » r. In the case of WMS, the probe beam modulation is generally 
treated as instantaneous frequency change: 



co i it) = — (o) 0 t + P sin co m t) -o) 0 + AF cos co m t 
At 



(2) 



where the maximum frequency deviation from the carrier frequency a>o is AF = fico m . In this model, 
dispersion effects cannot be considered, and the light intensity transmitted through the sample is 
expressed by It = Io exp(-oc), where a is the absorption coefficient. 

On the assumption that a(co) « 1, h{co) ~ h{co)[\ - a(co)], the a{co) can be expanded into the Taylor 
series around the carrier frequency coq. 
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Here cos co m t, cos co m t, etc., have been expanded with trigonometric identities for convenience in 
picking off the power at any frequency of interest. Because the instantaneous frequency varies 
sinusoidally, the Taylor series of It(co) is more or less automatically the Fourier series, too. This series 
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will not converge well for some AF [16]. However, if the frequency deviation AF is sufficiently small 
compared with F, higher-order terms can be neglected, and the WMS signal that is coherently detected 
at the frequency na> m is proportional to the n-th derivative of the sample absorption. 

In the present study, the assumption that AF « r was not valid under some experimental conditions. 
Thus we used the frequency domain approach to calculate the WMS spectral line shape. In the real 
system, the DFB laser intensity is not uniform over the FM range, and simultaneous amplitude 
modulation (AM) of the electric field occurs such that: 

E(t) = E 0 [l + M ski(co m t + y/)]exp[ico 0 t + ij3 sin(co m t)] (4) 

where M and y/ denote the AM index and the phase difference between AM and FM. By expanding the 
AM term in term of exponentials, i.e., Msm(co m t + yf) = M/2[cxip{i(co m t + y/)} - Qxp{i(co m t + y/)}], and 
replacing the exponential frequency terms as a summation over Bessel functions, one can rewrite 
Equation (4) as: 

oo 

E(t) = E 0 exp(ia) 0 t) £r„ exp(ina)j) (5) 

n = — oo 

where: 

r n = Z a k J n- k a 0 =l,a ±l = exp(± iyi) (6) 

k=-\ ^ l 

Equation (1) is the special case of Equation (5) where M = 0. 

When the modulated electric field E(f) is passed through a medium exhibiting the field amplitude 
attenuation S(co) and the optical phase shift tp(co), the transmitted optical field is then: 

oo 

E T (t) = E 0 exp(ico 0 t) £ exp(- 5 n - i<p n )r n exp(ina)j) (7) 

n=— oo 

where S n = S(a>o + nco m ) and tp n = tp(a>o + nco m ). Because the intensity / equals csoEE*/2 in SI units, the 
absorption coefficient a is twice the electric field amplitude attenuation coefficient S for weak 
absorption. The transmitted intensity can be expressed as a sum of Fourier components at frequencies 
corresponding to multiples of co m [13]. In the case of 2co m detection, the resulting signal is: 

oo 

/ = 21 0 [Re(z)cos 9 - Im(z)sin Q\ Z = £ r n+ / n _ x [exp(- (S n+i + 8 n _ x )) + i(<p n+1 - (p n _ x )] (g) 

«=-oo 

where 6 refers to the phase difference between the transmitted intensity and the reference signal in the 
lock-in amplifier. The in-phase (0 or 71) component corresponds to absorption, and the quadrature (±7i/2) 
components correspond to dispersion of the incident light. In the absence of absorption and dispersion, 
Z has a nonzero value due to residual amplitude modulation (RAM) [14]. In the case of 2co m detection, 
the background bias signal is: 

'ram=^ 2 'oCOs(6> + 2^ + ;t) (9) 
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Since /ram linearly depends on 7 0 , low frequency noise in 7 0 , such as fringe noise from optics and 
drifting of laser power and detector sensitivity, within the noise equivalent bandwidth of the lock-in 
amplifier is converted into FM/WM signal background noise. 

WM spectra are calculated from Equation (8) assuming that the absorption peak shape is described 
by the Voigt function. The Voigt function was calculated following the algorithm described by 
Humlicek [18] and Schreier [19]. The Doppler half width at half maximum (HWHM) was calculated 
from the temperature and molecular weight. To calculate homogeneous HWHM values, the pressure 
broadening factor was taken from the HITRAN database [20]. The FM index fi was estimated from 

= AF/co m . To calculate the AM index M, the intensity dependence on the laser frequency was 
measured, separate from the WMS experiments. Table 1 shows the parameters used in the spectral 
simulation. For the modulation frequency below 750 MHz, the AM-FM phase difference <// has been 
found to have a value of ca. nn + nil [14]. For the present study, negative dependence of the laser 
intensity on the laser frequency suggests that yj ~ -nil. In the WMS condition, i.e., co m « r, the 
dispersion term is negligibly smaller than the absorption term, thus the detection phase 6 ~ 0 for the 
maximum signal intensity. In all calculations, the values of y/ and 6 were fixed to be -tt/2 and 0, 
respectively. The signal intensity and noise are defined as peak-to-baseline amplitude and 
peak-to-peak variation of 2f signal with standard deviation of baseline, as shown in Figure 2. 



Table 1. Parameters used in the spectral simulation. 



Laser 


Target 


CO," 


ri a 
1 Doppler 


r b 

1 air 


/"self* 




1 


C0 2 


6,359.967 


0.00593 


0.0734 


0.1009 


0.423 


2 


N 2 0 


6,591.437 


0.00614 


0.0774 


0.1009 


0.305 


3 


CH 4 


6,046.953 


0.00933 


0.0627 


0.0820 


0.288 



" In units of cm" 1 . * In units of cm'Vatm. c In units of 1/cm" 1 . 



4. Results and Discussion 

For the detection of CO2, we used the 3vi + V3 combination band of the R(16) line at the absorption 
line center of 6,359.97 cm" 1 , which is free from interference by the other atmospheric species and has a 

23 2 1 1 

line strength of 1.73 x 10" cm molecule" cm" , according to the HITRAN database [20]. Figure 2 
shows experimental and simulated spectra of 1% CO2 at a total pressure of 15 kPa, diluted by N 2 . The 
experimental spectra were acquired with a lock- in time constant of 3 ms and filter slope of 24 dB/oct, 
which leads a noise equivalent bandwidth of the low pass filter to be 26 Hz. The simulated spectrum 
reconstructed the experimental spectrum well. Figure 3 shows the total pressure dependence of signal 
intensity and background noise. In all experiments, the modulation depth AF were adjusted to obtain 
the maximum signal intensity. The concentration of C0 2 was maintained at 1% in N 2 . In the simulation, 
we assumed that background noise was linearly dependent on M , according Equation (9), and the 
simulated noise intensities were normalized to minimize the residual to experimental values. As can be 
seen, the simulated results closely agree with the experimental measurements. The signal intensity 
profile of Figure 3 shows two opposing effects of pressure: decreasing spectral line centre intensities 
due to pressure broadening and an increasing absorber number density with increasing pressure. At low 
pressures (Doppler regime), the signal intensity increases linearly with pressure, while at high pressures 
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(Lorentzian regime) the signal intensity become almost constant due to the canceling of two opposing 
terms with pressure. Figure 4 shows experimental and simulated profiles of the signal to noise ratio 
(S/N) as a function of total pressure. The optimum total pressure in the detection of CO2 was found to 
be around 7 kPa. In the WMS, optimization of the signal intensity occurs at 
AF/HWHM = 2.2 in the Voigt function to calculate the properties of the 2/ line shape [21]. 

Figure 2. Experimental (solid line) and simulated (dashed line) 2f wavelength modulation 
spectra of 1% CO2 at around 6,359.97 cm" and a total pressure of 15 kPa diluted by N 2 . 
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Figure 3. Total pressure dependence of 2f wavelength modulation signal intensities 
(WMSI; •, - - ) and background noise (BGN; o, - - ) for 1% C0 2 at 6,359.97 cm" 1 
diluted by N 2 . Circles and lines indicate experimental and simulated results, respectively. 
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Because the Voigt HWHM is increased by pressure broadening with increasing pressure, the 
frequency deviation AF must be increased in order to maximize the WMS signal with increasing 
pressure. Since the AM index M depends linearly on AF, i.e., M/cQmfi (= M/AF) is constant, and the 
residual amplitude modulation /ram is linearly dependent on M according to Equation (9), the low 
frequency noise in Iq such as the fringe noise from optics and drifting of laser power and detector 
sensitivity was enlarged by factor of M at higher pressure, as shown in Figure 3. 

Figure 4. Total pressure dependence of S/N for 1% CO2 at 6,359.97 cm" 1 diluted by N 2 . 
Closed and open squares indicate experimental and simulated results, respectively. 
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In the case of N2O detection, the 3v3 overtone band of the R(15) absorption line at 6,591.44 cm" 

23 2 1 1 

with a line strength of 2.33 x 10" cm molecule" cm" was monitored. Figure 5 shows the 
experimental and simulated spectra of 1% N2O at a total pressure of 6 kPa. As for the previous 
experiment of CO2, the acquired 2/WM spectrum of N2O shown in Figure 5 is in excellent agreement 
with the simulated spectrum. Figure 6 shows the total pressure dependence of signal intensity and 
background noise for the 1 % N2O diluted in N2. The signal intensity is saturated above 20 kPa. 
Figure 7 shows experimental and simulated profiles of S/N as a function of total pressure. The 
simulated results agreed with the experimental results. The optimum total pressure for the detection of 
N2O is estimated to be around 6 kPa. 

For the detection of CH4, we used the R(3) lines of the 2v3 overtone band around 6,046.95 cm" 1 with 

212 11 

a line strength of the order 10" cm molecule" cm" [20]. As listed in Table 2, there are three 
absorption peaks at around 6,046.95 cm" 1 . Figure 8 shows the experimental and simulated spectra 
of 0.01% CH4 at around 6,046.95 cm" 1 and total pressure of 10 kPa. The three absorption lines were 
considered in the simulation. In the experimental spectra, only one absorption peak was observed 
at 6,046.95 cm" 1 due to the overlapping of the three absorptions. Figure 9 shows the total pressure 
dependence of signal intensity and noise for the 0.01 % CH4 diluted in N2. Figure 10 shows the 
experimental and simulated profiles of S/N as a function of total pressure. The simulated spectrum 
using the three absorption lines reproduces the experimental results well. The optimum total pressure 
of CH 4 detection at the 6,046.95 cm" 1 absorption peak was found to be around 15 kPa. The larger 
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optimum pressure for CH 4 detection compared with those for CO2 and N 2 0 would be caused by the 
broader nature of the absorption peak consisting of three closely separated absorption lines. 

Figure 5. Experimental (solid line) and simulated (dashed line) 2f wavelength modulation 
spectra of 1% N2O at around 6,591.44 cm" 1 and a total pressure of 6 kPa diluted by N2. 
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Figure 6. Total pressure dependence of 2f wavelength modulation signal intensities 
(WMSI; •, - - ) and back ground noise (BGN; o, - - ) for 1% N 2 0 at 6,591.44 cm" 1 
diluted by N2. Circles and lines indicate experimental and simulated results, respectively. 
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Figure 7. Total pressure dependence of S/N for 1% N 2 0 at 6,591.44 cm diluted by N 2 . 
Closed and open squares indicate experimental and simulated results, respectively. 
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Table 2. Absorption lines of R(3) of the CH4 2v3 band at around 6,046.97 cm" 



Line 


Position (cm" 1 ) 


2 11 

Integrated absorption cross-section (cm molecule" cm" ) 


1 


6,046.942 


8.17 x 10" 22 


2 


6,046.953 


1.00 x 10" 21 


3 


6,046.965 


1.32 x 10" 21 



Figure 8. Experimental (solid line) and simulated (dashed line) 2f wavelength modulation 
spectra of 0.01% CH4 at around 6,046.95 cm" 1 and a total pressure of 15 kPa diluted by N2. 



3. 
>> 

CO 
CD 

H — > 

CO 




-0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 

Av 0 (cm -1 ) 



Sensors 2010, 10 



4696 



Figure 9. Total pressure dependence of 2f wavelength modulation signal intensities 
(WMSI; •, - - ) and background noise (BGN; o, - - ) for 0.01% CH 4 at 6,046.95 cm" 1 
diluted by N 2 . Circles and lines indicate experimental and simulated results, respectively. 
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Figure 10. Total pressure dependence of S/N for 0.01% CH 4 at 6,046.97 cm" 1 diluted by N 2 . 
Closed and open squares indicate experimental and simulated results, respectively. 
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To estimate the limit of detection (LOD) of the present system, we changed the mixing ratio of 
sample gases, while keeping the optimum value of the total pressure fixed. The measurements were 
performed with a noise equivalent bandwidth of 26 Hz. From a fit of the signal intensity as a function 
of CO2 mixing ratio at 7 kPa, we found, for a signal to noise ratio of 2, the LOD of [CC>2] 7kPa = 24 ± 2 ppmv. 
In the case of N2O, the LOD at 6 kPa was estimated to be 7.6 ± 0.4 ppmv. For CH4, we estimated the 
detection limit and NEAS using the R(3) peaks at 15 kPa to be 60 ± 4 ppbv. The LODs for the three 
gases had an almost linear dependence on the absorption cross sections. The line strengths of CO2, N2O 
and CH4 of the used transitions vary slightly with temperature (0.23, 0.25 and 0.39%, respectively, by 
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temperature change of 1 K around room temperature), however, the LODs are unaffected by 
temperature change. Parkes et al. demonstrated the detection of N 2 0 by diode laser cavity ring-down 
spectroscopy using the same 3v3 overtone band [22]. The limiting sensitivity was 23 ppmv at 1 atm 
with mirror reflectivities of R = 0.9998, which is comparable with our 2/WMS value of 7.6 ppmv. The 
noise-equivalent absorption sensitivity (NEAS) per scan for a signal to noise ratio of 1 was estimated 
using following Equation (10): 

NEAS = ^^min 

iJnebw/n ( 10 ) 

where Afe m in, I, NEBW and N denote the absorbance at S/N = 1, the effective optical length of 29.91 m, 
the noise equivalent bandwidth of 26 Hz and the average number of 20 times, respectively. In the 
calculation of Abs m \ n , the absorption line shapes were assumed to be described as the Voigt functions. 
The estimated values of NEAS for C0 2 , N 2 0 and CH4 were found to be (1.3 ± 0.1) x 10" 8 , (4.8 ± 0.3) x 10" 9 , 
and (5.6 ± 0.4) x 10" 9 cm" 1 Hz" 1/2 , respectively. The higher value of NEAS for CO2 would be caused by 
the larger value of M/cOmfi, as shown in Table 1, which results the larger effect of the residual 
amplitude modulation. 

5. Conclusions 

We have described the highly sensitive detection of CO2, N2O and CH4 by 2/WMS, using tunable 
near-infrared diode lasers in the near- infrared region. The observed 2f WM spectra were analyzed as a 
function of the total pressure. At the optimum total pressure for the detection of each gas, as 
determined by experimental and theoretical approaches, the limits of detection in the present system 
were determined. Detection of atmospheric N2O (320 ppbv) at the 3v3 overtone band is not sufficient 
with our current WMS spectrometer. However, the necessary sensitivity could be delivered by reducing 
the optical fringe noise and the drifting of laser power using external modulators instead of direct laser 
modulation. 
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